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' Previously defined spin-dependent quark densities that are matrix elements of 

' specific density operators in proton states of definite spin-polarization generally have 

an infinite variety of non-spherical shapes. The present application is concerned with 
both charge and matter densities. We show that the Gross & Agbakpe model nucleon 
harbors an interesting variety of non-spherical shapes. 
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g ■ I. INTRODUCTION 

m , 

O . Recent data [l|, |^ showing that the ratio of the proton's electric and magnetic form factor 

^ \ Ge/Gm, falls with increasing momentum transfer for 1< < 6 GeV^ and equivalently 
that QF2{Q'^) / FilQ"^) is a approximately constant have created considerable attention. This 
behavior indicates that the sum of the orbital angular momentum of the quarks in the proton 
is non- vanishing 0, 0, S • 
^1 It is natural to consider the connection between orbital angular momentum f/) and the 

^ I shape of the nucleon, and one of us showed 0], using the proton model of Ref. that the 
rest-frame ground-state matrix elements of spin- dependent density operators reveal a host 
^ . of non-spherical shapes. The use of the spin-dependent density operator is the key feature 
^ I that allows the detailed connection between orbital, spin and total angular momentum 
to be revealed in quantum systems. Matrix elements of the non-relativistic spin-density 
operator have been measured in condensed matter systemsj^, revealing the orbital angular 
momentum content of electron orbitals. 

In the model of the relativistic nature of the quarks enters via their lower components of 
Dirac spinors that are part of the wave function. It is natural to ask if is it necessary that the 
presence of relativistically moving quarks always causes a non-spherical shape (as determined 
by matrix elements of spin-dependent density operators). We shall not address this general 
question here, and consider only the model dependence of the shape. In particular. Gross & 
Agbakpe (GA)|i3 have claimed to construct a manifestly covariant nucleon wave function 
that contains only / = such that "t/ie nucleon would still be spherical" and also describe 
the measured nucleon electromagnetic form factors. We show that this contention is not 
correct by evaluating the matrix elements of the spin- dependent density operator using the 
GA wave function. 
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II. SPIN-DEPENDENT DENSITY OPERATORS 

We begin by explaining how the shapes of a nucleon are exhibited by studying the rest- 
frame ground-state matrix elements of spin- dependent density operators. The usual charge 
density operator in non-relativistic quantum mechanics is given by 

P(r) = E 7^(r - r.) (1) 

where Cj/e is the charge of the i'th particle (in units of the proton charge) and its position 
operator. Matrix elements of this operator yield the charge density of a system. Suppose 
the particles also have spin 1/2. Then one can measure the probability that particle is at a 
given position r and has a spin in an arbitrary, fixed direction specified by a unit vector n. 
The spin projection operator is (1 -|- cr ■ n)/2, so the spin-dependent density operator is 

p(r,n) = ^^5(r-r,)^(l + ^-n). (2) 
. e z 

The spin-dependent density allows the presence of the orbital angular momentum to be 
revealed in the shape of the computed density. Matrix elements of the spin-density operator 
of Eq. have been measured in condensed matter systemsj^, revealing the orbital angular 
momentum content of electron orbitals. 

It is worthwhile to consider a simple example of a single charged particle moving in a 
fixed rotationally invariant potential in a state of quantum numbers: n, 1, 1/2, m we find 

-R^ (r) 

p(r, n) = (n, 1, 1/2, m |p(r, n)| n, 1, 1/2, m) = — !liM/^ (m|l + 2cr ■ r n ■ r — cr ■ n|m). (3) 

Suppose n is either parallel or anti-parallel to the direction of the proton angular momentum 
defined by the vector s. The direction of the vector s defines an axis (the "z-axis"), and the 
direction of vectors can be represented in terms of this axis: s ■ r = cos 9. With this notation 
p(r, n = s) = Rn^i i/2{r) cos^ 9, p(r, n = — s) = -R^ 1 1/2('^) sin^ 9 and the non-spherical shape 
is exhibited. The average of these two cases is a spherical shape as is the average over the 
direction of s. The spherical shape claimed in (lO| arises from taking an average over the 
direction of n. Averaging over this direction necessarily buries the effects of orbital angular 
momentum responsible for non-spherical shapes. 

The density we have discussed so far is defined in terms of position, but when quantum 
field theory applies it is convenient to define similar operators that give the probability for a 
particle to have a given momentum, K, and a given direction of spin, n. The field-theoretic 
version of the spin-dependent charge density operator is 

po(K, n)=J ^/^■^^l,{v)0{^' + 7 ■ n75)^(0), (4) 

where O is Q/e, the quark charge operator in units of the proton charge. Here we shall also 
consider the case in which O = 1. Its matrix element along with Eq. gives the spin- 
dependent matter densities. The quark field operators are evaluated at equal time. There 
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are no gluon fields in the relativistic constituent quark models used here. The question 
of color gauge invariance of QCD will be taken up elsewhere. The matrix element of this 
density operator in a nucleon state of definite total angular momentum defined by the unit 
vector s, |\l/s) is 

po(K,n,s) = (vl>,|pc,(K,n)|vl>,), (5) 

where the subscript O = Q, 1 specifies the operator used in Eq. 

The most general shape of the proton, obtained if parity and rotational invariance are 
upheld is 

Po(K, n, s) = Aa{K^) + Ba{K^)n ■ s + Co{K^)n ■ K s ■ K, (6) 

with the last term generating the non-spherical shape. Any wave function that yields a 
non-zero value of the coefficient Co(K^) represents a system of a non-spherical shape. 

In Ref. 0] the spin-dependent charge density of a relativistic three-quark constituent 
model was evaluated with the result (for the proton): 

p^(K,n,s) = p(Ji )^(l + n-s + 7(J0(l-n-s + 2K-nK-s)) (7) 

with 



p{K) = J d^k^'^{k,K){E{K) + m), 7(A') = 



E{K) - m 
E{K) + m' 



with K) the three-quark wave function and m the constituent quark mass specified in 
0]. The result for the neutron is 

p^(K,n,s) =p(ir)l(l-n-s + 7W(l + n-s-2K-nK-s)) (9) 

The shape for a given value of K is determined by the ratio 7 (-ft') which reaches j3| a 
value of 0.6 for K = 1 GeV/c. This implies considerable non-sphericity. The probability 
that a given value of K is determined by the function K'^p{K), displayed in Fig. 1 of [2|. 
The most likely value of -fC ~ 0.25 GeV/c corresponding to jlK) = 0.16. 

Some special cases of Eq. ((Tj) are interesting. Suppose the quark spin is parallel to the 
proton spin, n = s, then pQ{K,n = s) = p{K) (1 -|- 7(i^') cos^ 6'). For small K the shape 
is nearly spherical, but for large K the cos^ 6 term becomes prominent. On the other 
hand, the quark spin could be anti-parallel to the proton spin, n = — s. Then we find: 
p{K, n = — s) = p{K)'~f{K) sin^ 6, and the shape is that of a torus. We may also take 
the quark spin perpendicular to the proton spin n ■ s = 0, so that pQ(K,n -8 = 0) = 
p{K){l + 'y{K))/2 + 7(i^) sin6'cos^(cos0nx + sin^n^), to display the dependence on the 
azimuthal angle. In each case, the non-spherical nature arises from the term proportional 
to 7 (-ft') caused by the lower components of the Dirac spinor that is part of the proton wave 
function. For the neutron one obtains a doughnut shape for n = s, and a peanut for n = — s. 
Plots of these shapes have appeared in Ref. 01 as well as in the NY Times 
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FIG. 1: Matter distributions. The spin direction s is taken as vertical, z, with n = s. Right column 
model of 0] , left column model of First row K = 250 MeV/c, second row K=l GeV/c, third 
row K — > oo. 

For the matter distribution (O = 1) we find 

Pi(K, n, s) = p(K)-(3 + n ■ s + 7(K)(3 - n ■ s + 2K ■ nK ■ s)). (10) 

For a fixed value of K, this density looks more like a sphere than that of Eq. ((Tj). For example, 
if n = — s one finds pi(K, n) = p{K) {l +7sin^6'^, which is mainly spherical instead of 

toroid. Similarly if n = s one finds pi(K,n) = p(i^') (2 + 7 cos^ ^). For the case that 
n- s = 0, with n in the y-direction we find pi(K, n) = p{K) [3/2 + 7(3/2 + cos ^ sin 6* cos 0)]. 



III. MODEL OF GROSS & AGBAKPE 

In Ref. explicit effects of quark spinors are avoided by removing a factor of the 
single-quark free propagator that would appear in a solution of the Bethe-Salpeter equation. 
The GA nucleon wave function consists of two terms which involve either a scalar or vector 
diquark. The diquark covariant polarization vector is rj^. The relativistic impulse approxi- 
mation is used to evaluate the matrix element of the electromagnetic current operator. In 
the evaluation of the elecromagnetic matrix element, JT/" for a nucleon of isospin / and 
four-momentum transfer q is three times the integrated current of a single off-shell quark. 
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computed with the spectator diquark system on mass-sheU: 

1 



where P± = P ± |g. The CQ current, jj, is specified in [lO|. In deriving Eq. (fTT|) the 
polarization vectors of the diquark are summed over, using Yl,r]1^*1^' = — ~9^^' + 
.The quantity p is the on-shell diquark four-momentum (p^ = m^), and u(P, s) is the 
spinor of the nucleon. The scalar functions '?/'o,i describing the nucleon with scalar (0) and 
vector (1) diquark systems are chosen to depend only on the variable x = ((^ ~ ""^s)^ — 
(P — pY) / {Mms). For completeness we present the functions ilJoy. 

No 

= ™.(A-2 + x)(ft-2 + x)' 

where (5i and P2 are range parameters in units of Mrris, Nq is a normalization constant and 



We are concerned only with rest-frame matrix elements. In this frame x = '^{Es{i^) — 1), 



where Es{k) = a/T+k^ where k = y 1 + p^/m^. 

GA claim that their wave function depends only on the magnitude of p, and therefore 
is spherically symmetric. This is not correct. While the functions ipo^i depend only on the 
magnitude of p, the nucleon wave function contains the polarization vector r]'^ and this yields 
a non-spherical shape as defined by taking the matrix elements of the operator pc)(K,n). 
Before proceeding we note that the relation between the di-quark momentum p of GA is 
simply the negative of the quark momentum K of p = — K. 

The expectation value of the operator po{^, n) in the nucleon of total angular momentum 
s of the GA model (corresponding to Eq. (0)) is defined as pQ^(K,n,s). Once again the 
result is three times the integrated current of a single off-shell quark, calculated with the 
spectator "diquark" system on mass-shell. A straightforward evaluation yields 

pg^(K,n,s) = --i-^c7[0o(7o + n-775)|^o(K)2| 



+ 9^.^i^.7;.(7" + n ■ 775)7. - g'^j \M^)m (14) 

where U is the Dirac spinor for a nucleon of total angular momentum s at rest, and i?(K^) = 
JtC^ + ml. For the charge density: Oq = 2/3(proton), -1/3 (neutron), TjOiTj=0 (proton),-l 







FIG. 2: Matter distributions. The spin direction s is taken as vertical, z, with n 
column model of Q , left column model of [lo| . 
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(neutron). For the matter density: Oq = 1, tjOiTj = 3. Thus we explicitly find 

1 



GAp, 



MAn, 



pf^(K,n,s) 



K,n,s) 



K,n,s) 



2E(K2) 




There are a number of noteworthy features. The non-spherical nature of the distributions 
arise only from the vector di-quark part of the wave function. This does not contribute to the 
proton charge and spin-dependent charge densities, leading to a spherical (but polarization 
dependent) result. The charge distribution of the neutron and the matter distribution of 
the proton and neutron are each non-spherical with the degree of non-sphericity controlled 
by the terms TZ and 

Note that Eq. (fT3j) corresponds to a charge distribution of the proton because the charge 
operator is of the explicit form: 1/6 + T3/2 in That work accounts for the effects of 

the pion cloud 0, 01 in an approximate manner. In general, charge would be carried by 
both the quarks and the pion. In that case the quark charge operator would take the form 
/+ + /„r3, and the vector diquark term would contribute a non-spherical density proportional 
to 3/+ — /_. In general, quark hues are dressed by loops and the dressing can be applied to 
the quark charge density as shown in ET 
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GA 
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FIG. 3: Matter distributions. The spin direction s is taken as vertical, z, with n • s = 0. The 
direction of n is in the page {y direction). Right column model of 0], left column model of [lo| . 
First row K = 250 MeV/c, second row K=l GeV/c. 

IV. COMPARING SHAPES 

We turn to the numerical evaluation and display of these shapes, focusing on a comparison 
of the matter distributions of the Miller and GA models. Specific numerical evaluations 
of the shapes requires knowing the value of nig, which is not specified in Ref. Here we 
note that in the non-relativistic limit 7^ ~ 1 and the quantities K^/m^ and 7 play the same 
roles. In that case, we simply set K^/m^ = 7. We also consider a case when the quarks are 
moving relativistically with K = 1 GeV/c. Then 7 = 0.6. The relation between K and 7 
is K = 27727/(7^ — 27+1). It is reasonable to associate the values of 2m with the diquark 
mass nis. In this case, the value 7 = 0.6 corresponds to K"^ /ml = 3.75, x = 2.36, TZ = 0.53. 
These numerical values are used to obtain the figures shown below. 

The first situation we consider is the one in which the quark spin is parallel to the nucleon 
spin (total angular momentum). In the model 0] t he q uark probability is larger at the top 
than at the side by a factor of 1 + 7. In the model [lO| the corresponding factor is 1.24 for 
quarks of momenta 1 GeV/c. See Fig. 1. In the limit of A' ^ cxd the ratio would be 2 in 
both models. Thus the peanut configuration is hiding within the nucleon of , albeit with 
very small probability. 

Next we take the quark spin anti-parallel to the nucleon spin: n = — s, see Fig. 2. In this 
case the non- sphericity of the two models is very similar. Both have a toroid shape, with 
the model having more pronounced effects. 

The final situation we consider is that in which the quark spin is perpendicalar to that 
of the nucleon spin, see Fig. 3. We take n ■ s = with n pointing along the right side of the 
figure. Once again each model nucleon has a significant deformation. 
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V. SUMMARY 

The nucleoli is far from round in each of the models considered, and this arises from the 
relativistic nature of each. The deviation from a spherical shape is associated physically with 
motion of spin 1/2 quarks moving relativisitically within the nucleon, and mathematically 
with the non- vanishing of the term Cc)(K^) of Eq. The general nature of Co(K^) is a 
subject for future work. 
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